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ABSTRACT 

We present a theoretical analysis of galaxy-galaxy lensing in the context of halo models 
with CDM motivated dark matter profiles. The model enables us to separate between 
the central galactic and noncentral group/cluster contributions. We apply the model 
to the recent SDSS measurements with known redshifts and luminosities of the lenses. 
This allows one to accurately model the mass distribution of a local galaxy population 
around and above L*. We find that virial mass of galaxy is M200 = (5 — 10) x 
10 h Mq depending on the color of the galaxy. This value varies significantly with 
galaxy morphology with for late types being a factor of 10 lower in v! , 7 in g' and 
a factor of 2.5-3 lower in r', i' and z' relative to early types. Fraction of noncentral 
galaxies in groups and clusters is estimated to be below 10% for late types and around 
30% for early types. Using the luminosity dependence of the signal we find that for 
early types the virial halo mass M scales with luminosity as M oc T, 1A±0 - 2 in red bands 
above L+. This shows that the virial mass to light ratio is increasing with luminosity 
for galaxies above L*, as predicted by theoretical models. The virial mass to light 
ratio in i' band is 17(45)/iMq/Lq at for late (early) types. Combining this result 
with cosmological baryon fraction one finds that 70(25)%/i~ 1 Yif2 m /12f2f, of baryons 
within r2oo are converted to stars at L+, where Tj is the stellar mass to light ratio in i' 
band. This indicates that both for early and late type galaxies around a significant 
fraction of all the baryons in the halo is transformed into stars. 



Weak lensing by matter along the line of sight between the 
source and the observer shears the images of the background 
galaxies, inducing ellipticity distortions (see Bartelmann & 
Schneider 2001 for a review of weak lensing). Although away 
from rich clusters the effect is too small to be detectable 
for individual galaxy lenses, it can be measured statistically 
as a function of relative separation from the galaxy. This 
requires averaging over the tangential ellipticities of all the 
background galaxies relative to the lens and over all the 
lenses (galaxies which are in the foreground). Until recently 
this averaging, named galaxy-galaxy (g-g) lensing, was done 
as a function of apparent angular position in the sky (Tyson 
et al. 1984, Brainerd, Blandford, & Smail 1996, Hudson et al. 
1998, Fischer et al. 2000, Smith et al. 2001), so a signal at a 
given angular separation could be either coming from a small 
radial distance of a nearby lens or from a large distance 
of a far lens. This made the theoretical interpretation of 
the data rather involved. First attempt to use distances was 
by Wilson et al. (2001), which however only had limited 
photometric information and so could only obtain reliable 
distances to early type galaxies. 

Recent study of galaxy-galaxy lensing by SDSS collab- 
oration (McKay et al. 2001) is a significant step forward in 
the study of galaxy-galaxy lensing. The spectroscopic sam- 



ple of more than 35,000 lensing galaxies and 3.6 million back- 
ground galaxies is large enough to allow one a detailed study 
of the relation between mass and light for several luminosity 
bands and morphological types. Since the distances for lens 
galaxies are known one can study the strength of the sig- 
nal as a function of proper radial separation from the lens. 
In addition, because the survey is shallow, the redshift dis- 
tribution of background galaxies is known from the deeper 
spectroscopic surveys (e.g. Lilly et al. 1995). In combina- 
tion with above this means that the mean critical density 
is known for every lens, so one can average over the proper 
projected mass density rather than the shear itself. This fact 
greatly simplifies the theoretical analysis, since one can now 
measure the actual projected density as a function of galaxy 
luminosity and proper radial transverse distance from the 
galaxy. In addition, since the lens sample is at low redshift 
(mean z ~ 0.1) redshift evolution is small and k-corrections 
are relatively reliable in red bands, further simplifying the 
theoretical interpretation. 

In this paper we want to connect the SDSS observa- 
tional results to the theoretical models in the context of our 
current understanding of galaxy formation models within 
the CDM paradigm. We will model the dark matter ha- 
los with CDM type of halo profiles (NFW profile; Navarro, 
Frenk, & White 1997), where the slope is gradually changing 
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from the inner slope between -1 to -1.5 to the outer slope of 
-3. Although other profiles have been proposed that differ 
significantly from NFW in the inner parts of the halo, they 
agree well with NFW in the outer parts (Klypin et al. 2001). 
NFW profile should be contrasted to the truncated singular 
isothermal sphere with a constant slope -2 out to a fixed 
radius which was often used used in the past work on g-g 
lensing. Our main goal is to determine the virial mass of the 
halo and its relation to the luminosity of the galaxy. This is 
important for theoretical models of galaxy formation, since 
it is usually assumed that only baryons within the virial ra- 
dius are able to condense and form stars. By determining 
the virial mass for a given galaxy luminosity one can thus 
directly determine the efficiency of star formation in typical 
galactic halos. Luminosity dependence of the galaxy-galaxy 
lensing signal allows one to determine this as a function of 
halo mass, while morphology subsamples can determine it 
as a function of morphological type. By comparing the ra- 
dial dependence of the signal with the theoretical models 
one can also determine the fraction of galaxies that live in 
larger halos such as groups and clusters. This is another pa- 
rameter that can distinguish between the different galaxy 
formation models. Finally, comparison of high density sam- 
ple to the field sample allows one to study the effects of 
dense enviroments, such as tidal stripping, on the dark mat- 
ter profile and mass to light ratios. Thus g-g lensing allows 
one to make detailed tests of the galaxy formation models 
(e.g. Kauffmann et al. 1999, Benson et al. 2000, Somerville 
& Primack 1999). In principle g-g lensing can also allow a di- 
rect determination of the dark matter halo profile, although 
as we will show here at present the data do not have enough 
power to strongly constrain this. 

An important issue in the theoretical analysis of g-g 
lensing is how to separate the contribution from the indi- 
vidual halo of the galaxy from that of the neighbouring 
galaxies or larger mass concentrations such as groups and 
clusters. The former should dominate on small scales while 
the latter on large scales. Previous work used the additional 
information obtained from the galaxy clustering to remove 
this contribution, but in general such analysis relies on the 
assumption that all the mass is associated with galaxies, 
which is invalid on scales where groups and clusters become 
dominant. Most of the mass in these systems is in a diffuse 
form and only about 5-10% is expected to be attached to 
individual galaxies (e.g. Springel et al. 2001, Ghigna et al. 
2000, Bullock et al. 2001). This means that these systems 
cannot be modeled using the galaxies as the mass tracers, 
which only accounts for a small fraction of the total mass. In 
fact, theoretical models presented in this paper suggest the 
contribution from these group and cluster halos dominates 
the signal on scales above 200/i~ 1 kpc (see also Seljak 2000) 
and has to be carefully modeled to account for it. 

Given that corrections to the profile on large scales are 
difficult to determine from the data directly one has to turn 
to theory for guidance. Realistic theoretical modeling of g- 
g lensing must combine galaxy formation models and dark 
matter models. One way is to use semi-analytic or hydrody- 
namic models of galaxy formation and combine them with 
N-body simulations (Guzik & Seljak 2001; White, Hernquist, 
& Springel 2001). This has the advantage of having a real- 
istic distribution of dark matter and galaxies, but suffers 
from the lack of force and mass resolution. In addition, this 



approach by itself does not allow for a fast exploration of 
parameter space. The alternative approach is to use the 
recently popularized halo model applied to galaxies (Sel- 
jak 2000, Peacock & Smith 2000, Scoccimarro et al. 2001, 
Berlind & Weinberg 2001), which takes into account both 
the individual halo profiles (for galaxies either at the cen- 
ters of the galactic halos or distributed within larger groups 
and clusters with a specified radial distribution) and corre- 
lations between the galaxies. This approach applied to the 
galaxy-dark matter correlations (which fully determines the 
g-g lensing signal) was shown to give the same results as 
the simulations in the regime of applicability (Seljak 2000). 
However, halo model is analytical, does not suffer from the 
resolution issues and provides a more physical interpretation 
of the results. Thus one can parametrize the model with the 
quantities one wishes to extract from the data and deter- 
mine these directly. It also allows for a rapid exploration of 
the parameter space without the need to rerun cosmological 
simulations or to repopulate the halos with galaxies using 
semi-analytic galaxy formation models. In this sense the halo 
model provides a natural link between the observations and 
the theoretical models of galaxy formation. 

Even though we will use the halo model to analyze g- 
g lensing the main features can be understood without it. 
Particularly robust conclusions are possible using the low 
density sample, where the clustering and group/cluster con- 
tribution can be neglected. In this case the signal can be 
taken simply as a projected radial mass profile of the halo 
averaged over the mass distribution of the halos determined 
by the galaxy sample. For the full and high density sam- 
ple the contribution from groups and clusters can no longer 
be neglected, as it dominates the signal on scales above 
200/i _1 kpc. This allows one to determine the fraction of 
galaxies in groups and clusters as well. 

The outline of this paper is as follows: in section §2 
we review the basic theory and halo model as applied to 
g-g lensing. In §3 we discuss the influence of the parame- 
ters introduced on the observed g-g lensing signal. We fo- 
cus on the relative contributions to g-g lensing from galaxy, 
group/cluster and clustering terms and explore which pa- 
rameters the observations are most sensitive to. In §4 we 
apply the model to the data to determine several of the 
model parameters. Interpretation of the results and conclu- 
sions are presented in §5. 



2 CONTRIBUTIONS TO G-G LENSING FROM 
THE HALO MODEL 

In the halo model it is assumed the galaxies form in collapsed 
dark matter halos. The halos are clustered among them- 
selves and have a specified density profile. The contribution 
to g-g lensing comes either from the dark matter profile of 
the halo the galaxy is sitting in (one-halo or Poisson term) 
or from the clustering of all the other halos around a given 
galaxy (halo-halo term). As we show below the second con- 
tribution is small compared to the first one on scales below 
l/i _1 Mpc of interest here. The one-halo term has two con- 
tributions. First is from the dark matter around the galaxy 
itself that depends on the dark matter halo profile. Second is 
from the dark matter in groups and clusters the galaxy may 
belong to, but is not at their center. In the latter case one 
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must also specify the radial distribution of galaxies inside 
these halos. The simplest case is to assume it is the same as 
the radial distribution of the dark matter (see e.g. Hoekstra 
et al. 2001 for some justification in the case of observations 
and Ghigna et al. 2000, Springel et al. 2001 for the case of 
simulations) . 

The relative contribution from the two one-halo compo- 
nents depends on the fraction of galaxies residing in groups 
and clusters. In general this fraction varies with luminos- 
ity, type and other selection criteria. Particularly simple is 
the low density sample in SDSS data, which consists mostly 
of the field galaxies and where the group/cluster fraction is 
likely to be small. When this is not the C£LS6, £LS for example 
for the high density selected sample or early type sample, 
one must also address the question of whether the halo pro- 
files around galaxies that are within larger halos differ from 
the profile of equal luminosity galaxies in the field. Even if 
the galaxies that end up in groups and clusters do not have 
different star formation efficiencies, so that their luminosity 
is a good indication of the halo mass, they may still loose 
some dark matter after merging. This depends on the amout 
of tidal stripping, which itself depends on the radial orbit of 
the galaxy. Simulations indicate that the tidal radius should 
be defined at the orbit pericenter, which is the point of the 
closest approach to the center of the cluster (Ghigna et al. 
2000). Since the orbits are very eccentric (typical ratios 1:5) 
the pericenter can be significantly smaller than the instanta- 
neous galaxy position. Still, these simulations indicate that 
the halo profiles of most of the halos within halos remain un- 
affected out to where the density of the larger halo matches 
that of the smaller halo (Ghigna et al. 2000). Since it is the 
radial halo profile that is relevant for g-g lensing we will 
model the galaxies within clusters as a superposition of a 
smooth cluster profile (weighted by the radial distribution 
of galaxies) and an individual galactic halo profile, which is 
assumed to be the same as for the equal luminosity galaxy in 
the field. We discuss further below how this assumption can 
be tested with g-g lensing itself. We do not need to model the 
transition from the subhalo to the cluster, where there may 
be departures from this simple model, since one is observing 
a projected density with g-g lensing. Given that clusters are 
significantly larger than galactic halos the cluster dominates 
the signal at radii smaller than the 3-d transition point. 



E crit7i (7?) = AE(iJ) = E(fl) - £(#). 



(2) 



Here R is the radial distance from the galaxy and E(i?) is 
the mean surface density within R. AE(i?) is the quantity 
that is observable, so all of our predictions will be compared 
to it. Critical surface density depends on the distances to 
the lens and source galaxies, 



Ecrit — 



D s 



A-kG DlDls' 



(3) 



where Dl and Ds are the angular diameter distances to 
the lens and source, respectively, and Dls is the angular 
diameter distance between the two. In the case of SDSS wc 
know the redshift of the lens, so Dl is known, while Ds 
and Dls can only be determined in an average sense. Since 
SDSS survey is so shallow, the redshift distribution of back- 
ground galaxies is well determined from other deeper red- 
shift surveys such as CNOC (Lilly et al. 1995). This means 
that there is little error associated with the redshift distri- 
bution of background population, which is not the case for 
deeper lensing surveys. While there is some uncertainty in 
the transformation from redshifts to distances this is gen- 
erally a small effect given the low redshifts of background 
population (z ~ 0.3 — 0.5). Here we will ignore it, assum- 
ing distances using cosmological constant model with matter 
density Q m = 0.3. 

To relate the observations to theory we interpret the 
above expressions in an average sense by averaging over all 
the galaxies and their surrounding dark matter. In this case 
the surface density E is related to the galaxy-dark matter 
correlation function (Guzik & Seljak 2001) 



E(i?) 



4dm[(H a + X 2 ) 1/2 )dx 



(4) 



where we dropped the unobservable constant term. Since 
we are interested in the small scales where galaxies do not 
evolve much we will express everything in proper coordi- 
nates. The low mean redshift of lens population (z ~ 0.1) 
and the fact that most of the galactic halos have formed at 
higher redshift means that the redshift evolution effects will 
be small if expressed in proper coordinates and relative to 
the matter density today We are implicitly assuming that 
most of the signal is arising near the galaxy, so we ignore 
the changes in the focusing strength along the line of sight. 



2.1 Theory 

We are interested in the distribution of the dark matter 
around the galaxies of a selected type, averaged over all the 
galaxies in the sample. We can quantify this as an excess of 
the dark matter density above the average as a function of 
the radial separation r. This is described by the galaxy-dark 
matter cross-correlation function, 

Cg,dm(r) = (dg(x)Sdm(x + r)) . (1) 

Here <5 g and <5d m are the overdensities of galaxies and dark 
matter, respectively 

Observationally the quantity we are interested in is tan- 
gential shear y t , which describes elongation of images per- 
pendicularly to the line connecting the image and the lens. 
This is related to the convergence k = E/E cr it, where E is 
the projected surface density in units of critical density, 



2.2 Halo model 

For a given correlation function £ gl dm we can use equations 
(2-4) to calculate AS. In the halo model it is easier to calcu- 
late the power spectrum P g ,dm(fc), which is just the Fourier 
transform of the correlation function. The halo model for 
power spectrum assumes the matter is in a form of isolated 
halos with a well defined mass M and halo profile p(r, M). 
The latter is defined to be an average over all halos of a given 
mass and does not assume all halos have the same profile. 
The mass is determined by the total mass within the virial 
radius r v i r , defined to be the radius where the mean density 
within it is <5 v i r times the critical density of the universe. 
Here we will use the value <5 v ir = 200 relative to the criti- 
cal density today. Often one uses the value defined by the 
spherical collapse model, which gives <5 v i r ~ 100 for ACDM 
model with fi m = 0.3. The latter gives about 15% higher 
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virial mass than the values presented here for a typical halo 
profile. 

We model the halo density profile in the form 



p[r) = 



(r/r s )-«(l +r/r a ) 3 + c 



(5) 



This model assumes that the profile shape is universal in 
units of scale radius r s , while its characteristic density p s at 
r s or concentration c = r v /r s may depend on the halo mass. 
The halo profile is assumed to scale as r~ 3 in the outer parts 
and as r° in the inner parts, with the transition between the 
two at r a . We fix the inner slope to a = — 1 (Navarro, Frenk, 
& White 1997), since SDSS g-g lensing is not sensitive to 
small scales given that the first bin is at R = 75/i _1 kpc. For 
this profile the typical concentration on galactic scales is 
expected to be C200 ~ 10 (Bullock et al. 2001, Eke, Navarro, 
& Steinmetz 2001). We explore the sensitivity of g-g lensing 
to this parameter below. 

The first contribution to the one-halo or Poisson term 
comes from the dark matter halos attached to the galaxy, 
where the galaxy is assumed to be at the halo center. If all 
the galaxies were of the same mass then this term would be 
just the Fourier transform of the dark matter profile itself. A 
generalization of this is to integrate over the luminosity dis- 
tribution assuming a relation between the galaxy luminosity 
L and halo mass M(L), 



peg 
g,dm 



(fc) = 



1 



dn 



(27r) 3 n J dL 



y[k,M(L)]dL, 



(6) 



where dn/dL is the known luminosity distribution of the lens 
sample and n is the total number of galaxies in the sample. 
The superscript eg stands for central galaxy. Here y(k, M) = 
p(k, M) I 'M , where p(k,M) is the Fourier transform of the 
radial halo profile of dark matter 



p(k,M) = / 4nr 2 drp(r,M)- 



i(kr) 



kr 



(7) 



Note that unlike in previous work on halo model we do not 
require the profile to stop at the virial radius, so there can 
be more mass associated with the halo than the virial mass 
itself. This is only possible if we do not require the mass 
function defined below to integrate to unity, since some of 
the mass can be part of the extended halo structure beyond 
the virial radius. This is the case for the mass function by 
Jenkins et al. (2001) that we use here. Typically we extend 
the halo to 2-3 times the virial radius and we verified that 
the exact cutoff position makes very little difference in the 
final results. As discussed above we assume the dark matter 
profile for a given luminosity does not depend on whether a 
galaxy is in a larger halo or not. In the extreme case that the 
galaxies within larger halos do not have any mass attached 
to them at all the virial masses of field galaxies should be 
underestimated by the fraction of galaxies in groups and 
clusters. Most of these galaxies are in fact in the field, so 
this fraction and the corresponding correction is small. 

Second contribution to the one-halo term comes from 
the dark matter around galaxies in groups and clusters, 
where the dark matter is not associated with the galaxy 
but with a larger halo. This term by definition includes only 
the galaxies that are not at the center of the halo they find 
themselves in, since such galaxies are already accounted for 
by the previous term. To obtain the g-g signal one must 



average over the mass distribution of halos and over the dis- 
tribution of galaxies within these halos. Since the halos come 
in a variety of sizes we must first introduce the halo mass 
function dn/dM, describing the number density of halos as 
a function of mass. It can be written as 

dn = £_ f( a ) dlnc7 ~ 1 ( 8 ) 
d\nM M JK ' d\nM ' V ' 

where p is the mean matter density of the universe, M is 
the virial mass of the halo and n(M) is the spatial number 
density of halos of a given mass M. We introduced a func- 
tion f(a), which has a universal form independent of the 
power spectrum, matter density, normalization or redshift if 
written as a function of rms variance of linear density field 
a. Jenkins et al. (2001) propose the following form (see also 
Sheth & Tormen 1999), 



f(a) = 0.315 exp[-| In a' 1 +0.61| 3 8 ], 



(9) 



which they argue is universal if mass used is within the ra- 
dius where overdensity is 200 of mean density (since we use 
a different definition of halo virial mass we must correct for 
the difference between the two). Note that this mass function 
does not account for all the mass density in the universe, so 
that some of the mass could be present in the outer parts of 
the halos. 

We must still specify how the galaxies populate the ha- 
los of different mass. For any galaxy sample this can be 
parametrized as the mean number of galaxies as a function 
of halo mass (N)(M) (Jing, Mo, & Boerner 1998). Here we 
work with (N) (M) averaged over the luminosity distribution 
of the sample, which is dominated by the galaxies around 
(see figure 1). We describe a model for (N)(M) below. 
The contribution from this term is 



g,dm 



(*) 



1 



(27r) 3 n 



f(a)dlna(N)y(k,M)y g (k,M), (10) 



where y g is defined as the Fourier transform of the radial 
distribution of galaxies and the superscript gc stands for 
group and cluster contribution. We will assume y g = y in 
most of this paper, but we also explore the sensitivity to this 
assumption in the next section. 

Finally, the halo-halo contribution to the power spec- 
trum describes the correlations between the lens galaxy and 
neighbouring halos. Here we assume it follows the linear 
power spectrum, except that halos can be biased relative to 
it, with bias b(a) a function of rms variance a. This assump- 
tion has been shown to give a good agreement with simula- 
tions of galaxy-dark matter correlations (Seljak 2000). Low 
mass halos are unbiased or mildly antibiased (b < 1, Sheth, 
Mo, & Tormen 2001). The power spectrum of this term can 
be calculated using (Seljak 2000) 



PgXv(k) = Plin(fc) 



/ 



I f f(a)d\na^b(a)y p g [k,M(a) 



f(a)dlnab(a)y[k,M(a)} 



(11) 



where p — for the field component and p = 1 for the 
group/cluster component. Here we implicitly assumed that 
(N)(M) includes also the central galaxy component, which 
means it will have a strong peak at the halo mass corre- 
sponding to the Li, luminosity (figure 1). Since we will show 
that this term is negligible on the scales below l/i _1 Mpc of 
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interest here the details of its modeling are not very impor- 
tant for our purpose. 

The total power spectrum is simply the sum of all con- 
tributions. Once this is computed we Fourier transform it 
to obtain the correlation function and integrate it along the 
line of sight (equation 4) to obtain the projected surface 
density E(i?). Another integration (equation 2) is needed to 
obtain the mean surface density within R and the observed 
quantity AS(_R). Since the redshift of lenses is low we do 
not include any redshift evolution effects. 

The above model left several functions unspecified, most 
notably the halo mass-luminosity relation and the galaxy 
occupation number (N)(M). We assume that there is a one 
to one correspondence between the halo mass and luminosity 
of the galaxy, which we model as a power law 



M _ (J^Y 



(12) 



where M 4 is the mass associated with L* galaxy. The slope 
P parametrizes the relation for a simple power law depen- 
dence. This is not the most general possibility and over a 
broad range of luminosity /3 is expected to vary, but since 
we are working over a narrow range of luminosities (a fac- 
tor of 10 at best) this parametrization should not be too 
bad and in any case, the data do not yet allow for a more 
general parametrization. Note that bright red galaxies have 
been removed from the sample, since it is likely that they 
do not fall on the same mass-luminosity relation. Both M* 
and (3 can be extracted from the data. We could also give 
some scatter to this relation, but in practice the scatter is 
dominated by the fact that we cannot choose very narrow 
luminosity bands for the analysis and so we will not pursue 
it further here. The distribution of halo masses has a width 
that depends on the luminosity distribution of the galaxies 
as shown in figure 1. 

The second unspecified function is the galaxy occupa- 
tion number (N)(M) for noncentral galaxies in groups and 
clusters. This is expected to grow with the halo mass for any 
galaxy population. Here we assume that the mean number 
of galaxies (of a given type) has a power law relation to the 
mass of the halo, 



(N)<xM% M > M cutoS . 



(13) 



We assume the power law relation is valid from some mini- 
mum halo mass M cuto ff , which should be at least a few times 
above the typical halo mass for a given luminosity sample, 
since by definition it is assumed these galaxies are not at 
the halo centers, where presumably another bright galaxy 
resides. We explore the sensitivity to the low mass cutoff 
A^cutoff and e below. We normalize the group/cluster con- 
tribution in terms of a fraction of galaxies a that reside in 
these larger halos. 

In figure 2 we show {N)/M versus M as calculated from 
semi-analytic simulations (Kauffmann et al. 1999) for two 
narrow intervals in luminosity, a brighter one around L* 
and a fainter one magnitude below it. One can see that in 
both cases the distribution is well modeled as a sum of a 
narrowly peaked distribution for the central galaxies and a 
linear power law component (N) oc M for the noncentral 
galaxies in groups and clusters. In our model we would pre- 
dict first contribution to be narrow because we have chosen 
a narrow bin in luminosity, which is strongly correlated with 



T3 




-22 -20 -18 
M (mag) 



Figure 1. Magnitude distribution of the SDSS sample in i' band 
(in units with h = 1) for the whole sample and the 4 luminosity 
subsamplcs separated by dotted lines. Also indicated is the L* 
magnitude (Blanton et al. 2001) as a vertical dashed line, which 
roughly corresponds to the mean luminosity of the whole sample 
((L) ~ L* = 2.05 X 10 10 /i~ 2 Lq). Note that even though the 
mean luminosity is close to L+ most of the lensing signal comes 
from the massive galaxies above L*, which explains the choice of 
luminosity subsamples. 



the halo mass for central galaxies. This is clearly seen to be 
the case in figure 2. Second contribution is a power law be- 
cause more massive halos also contain more subhalos of a 
given mass, which host noncentral galaxies. It is the latter 
contribution that enters in our model for (N)/M in equa- 
tion 6, while the first component is included by the galaxy 
luminosity function. Our model therefore qualitatively re- 
produces the main features of SAMs, but by characterizing 
the main ingredients with a few free parameters it allows us 
to be more general than any specific model. 

It is interesting to compare the fraction a of galaxies 
in groups/clusters for the two cases from figure 2. This is 
20% for the bright sample and 25% for the faint sample. 
This parameter can therefore vary as a function of lumi- 
nosity and is expected to decrease with L. This means that 
the group/cluster contamination is less important for bright 
galaxies, where it becomes more difficult to determine it ob- 
servationally because of the small number statistics. We will 
use this fact in the modeling below. 



3 DEPENDENCE OF G-G LENSING ON THE 
HALO MODEL PARAMETERS 

3.1 Halo mass and concentration 

We have introduced several free parameters in the model 
above and in this section we explore the dependence of the 
measured signal on them. The main parameter is the virial 
halo mass M* of a typical galaxy. The resulting g-g lens- 
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Figure 2. N(M) calculated from SAMs using GIF simulations 
(Kauffmann et al. 1999) using two narrow luminosity samples, 
-20.5 < I < -20 (dashed) and -22 < I < -21 (solid). Both the 
central galaxy component peaking at low halo masses and non- 
central group/cluster component are included here. The latter has 
N(M) approximately proportional to M, ie e = 1. 



Figure 3. Signal for a NFW halo profile with c = 10 of vary- 
ing mass. From bottom to top M = 10 n ,3 x 10 11 ,10 12 ,3 x 
10 12 , 1O 13 /i- 1 M . Also shown arc the observational data for the 
low density sample. 



ing signal for halos of masses between 10 11 - lO 13 fe _1 M 
is shown in figure 3. For reference we also show the SDSS 
data using the low density subsample, described in more de- 
tail below. Only the one-halo term is shown and we used 
NFW profile with c = 10 and a 5-function in the halo mass 
distribution. What is shown is therefore a simple projection 
of NFW profile, which has a well known analytic form for 
£(.R) (Bartelmann & Schneider 2001). The signal at a given 
scale is rapidly increasing with the mass of the halo, indi- 
cating that SDSS g-g lensing is a very sensitive probe of the 
virial mass of the halo. Masses around 10 12 h^ 1 Mq seem to 
be required to explain the field sample. 

The dependence on the concentration parameter c is 
shown in figure 4. We keep the virial mass of the halo con- 
stant, so the signal at large radii is similar for all the mod- 
els. At smaller radii the profiles start to differ and the more 
concentrated halos give a larger signal than the less concen- 
trated ones. Doubling the concentration index c changes the 
signal by 20-30% in the inner most bin at 75/i -1 kpc and less 
than that at larger radii. We will see below that since some 
contribution to the mass determination comes from the in- 
ner most bin there is some degeneracy between c and M, 
such that larger values for c give lower values for M. How- 
ever, this degeneracy is not complete and too steep profiles 
can be ruled out regardless of the halo mass. 

Since we cannot choose an infinitely narrow distribution 
in luminosity and have a high signal to noise a 5-function 
in the halo mass distribution is not very realistic. Instead 
one must integrate over the halo mass distribution to ob- 
tain the average profile. While the mean signal will corre- 
spond to the mean halo mass, the shape of the signal may 
be more affected, possibly complicating the determination 




200 
[h" 1 kpc] 



400 



Figure 4. Signal for a halo with NFW profile of fixed mass vary- 
ing concentration c from 3 to 16. 



of the concentration parameter. To investigate this possibil- 
ity we compute the signal for the luminosity distribution in 
figure 1, assuming halo mass is proportional to luminosity 
and comparing it to a single NFW profile corresponding to 
the mean mass of the sample. The result is shown in figure 
5. The two profiles are very similar, indicating that even for 
realistic luminosity distributions the shape of the halo can 
be reliably determined from the data if all the galaxies are 
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Figure 5. Comparison between a single NFW profile for M* 
(solid) and an average profile using the realistic luminosity distri- 
bution and assuming luminosity is proportional to mass (dashed). 

at the halo centers. These results show that if the galaxies 
are in the field then g-g lensing is a very robust probe of the 
halo virial mass and its profile. The field assumption should 
be valid for the low density sample discussed in more detail 
below. For the whole sample one must introduce also the 
contribution from groups and clusters, which is discussed 
next. 

3.2 Group and cluster contribution 

The group and cluster contribution to lensing is shown in 
figure 6 for e = 1. Here it is assumed all the galaxies are 
in groups and clusters (a = 1). It is evident that this con- 
tribution is significantly different from the central galaxy 
contribution. At small radii the group/cluster contribution 
is negligible, since by definition these galaxies are not at 
the halo center (the individual subhalo contribution from 
these noncentral galaxies is already included in the previ- 
ous term), the cluster component itself is very smooth and 
the galaxies are assumed to be distributed radially like dark 
matter and do not have a strong central peak. We remind 
the reader that for a uniform density the shear signal goes 
to zero by equation 2. The group/cluster contribution peaks 
around 200/i _1 kpc and then gradually decreases with radius, 
such that even at l/i _1 Mpc it is still 60% of the peak value. 
One can see that varying the low mass cutoff M cut off has 
some effect on the amplitude of the signal. For lower mass 
cutoffs one finds a lower average signal, because there are 
proportionally more galaxies in less massive halos. On the 
other hand, the radial dependence of this contribution is less 
affected. Since the radial dependence of this contribution is 
very different from the central galaxy contribution (figure 3) 
the two contributions can be robustly separated in the data 
even if the actual value of a has some uncertainty related 
to the low mass cutoff. Based on the results in figure 2 wc 



Figure 6. Group and cluster contribution to g-g lensing as a 
function of lower mass cutoff M cuto ff • Results shown here are for 
Mcutoff of 10 12 /i _1 Af Q (bottom curve), 3 X 10 12 /i _1 Mq (middle 
curve) and 10 1A h~ 1 Mq (top curve). We use c = 7 for galaxy and 
dark matter halos. 



choose M C u to ff to be 3 times the halo mass corresponding to 
the luminosity of the sample. For L* galaxies this is around 
2 x 10 12 /i _1 Mq. This cutoff is increased when we analyze 
brighter luminosity subsamples. 

In addition to the group low mass cutoff we also ex- 
plored the sensitivity of this component to the slope e of the 
galaxy occupation number (equation 13) and to the radial 
galaxy distribution y g (equation 6) . The results are shown in 
figure 7. Higher e gives more weight to more massive clusters 
and has a higher signal than a lower e. Similarly, more con- 
centrated galaxy distribution enhances the signal since now 
more galaxies are closer to the cluster center where the halo 
density is higher. Most observations and simulations predict 
that galaxy distribution does not deviate significantly from 
the dark matter, which has a mean concentration of around 
c = 5 — 10 for groups and clusters, depending on the low 
mass cutoff. Here we compare the cases with c = 3, 7, 10 
for the galaxy while using c = 7 for the dark matter. In 
all cases the overall amplitude changes somewhat depend- 
ing on the parameters, while the qualitative shape of the 
curve remains similar, peaking around 200/i _1 kpc and then 
decreasing with radius. This is distinguishable from the cen- 
tral galaxy component regardless of the specific choice of the 
parameters. In the following we assume e = 1 and y g = y 
with c — 7. 

To combine the central galaxy and group and cluster 
contributions we must assume parameter a, which is the 
fraction of galaxies in groups and clusters. Variation in this 
parameter is shown in figure 8. If the fraction a is small 
the signal on large scales is also small, whereas in the oppo- 
site case the large scale and small scale signals are com- 
parable. Also shown are the data from the whole SDSS 
sample, which show a relatively flat signal between 300- 
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Figure 7. Group and cluster contribution to g-g lensing as a func- 
tion of e=0.8, 1.0, 1.2 (short dashed, solid, long dashed, respec- 
tively) and using c = 7 for galaxy and dark matter. Also shown 
is variation of radial galaxy distribution parametrized with con- 
centration c = 3 (dot-short dashed) and c = 10 (dot-long dashed) 
for e = 1. 



lOOO/i - kpc. This cannot be explained by the central galaxy 
contribution alone, which drops below the observed points 
at R > 200/i _1 kpc. It also cannot be explained by all the 
galaxies being in groups and clusters, which would have an 
excess of signal at large radii. A contribution from noncen- 
tral galaxies in groups and clusters with a = 0.2 together 
with the central galaxy contribution provides a good de- 
scription of the data. 

The contribution from the correlations between the 
galaxies and the neighbouring halos (halo-halo term) is also 
shown in figure 8. This contribution is very small and can- 
not explain the strength of the signal on scales around 300- 
lOOO/i kpc. Instead the group and cluster contribution is 
required to explain the signal on these scales. Note that this 
does not mean that the correlations between the galaxies are 
not important, since the galaxies residing in groups and clus- 
ters are also correlated among themselves. It does however 
mean that the correlations between the galaxies in isolated 
halos are small and are instead dominated by the galaxies 
that reside in groups and clusters, whose contribution to 
g-g lensing is included in group and cluster term discussed 
above. The main difference between the two terms is that for 
the latter most of the mass in large halos is not associated 
with the halos of individual galaxies within them, which ac- 
count for only about 5-10% of the mass, but with a diffuse 
component of groups or clusters. In this respect our anal- 
ysis differs from the previous analysis of this data (McKay 
et al. 2001), where only the mass associated with the visible 
galaxies was accounted for in the correlation analysis. 



Figure 8. Central galaxy contribution for M* = 10 h~ Mq 
(solid), group/cluster contribution with a = 0.2 (short dashed) 
and a = 1 (long dashed) and halo-halo contribution to g-g lensing 
(dotted). Also shown is the SDSS data for the whole sample. 

4 APPLICATION TO THE SDSS DATA 

In this section we apply the formalism developed in previous 
sections to the SDSS data. We begin first by analyzing the 
luminosity dependence of g-g lensing, which provides us with 
a parametrization of the mass luminosity relation and with 
a typical mass M* of an L* galaxy. We proceed by analyzing 
the full sample as well as high and low density versions of it, 
which provide an alternative estimate of M* and the fraction 
a of group/cluster galaxies in the sample, as well as an esti- 
mate of mass loss in dense enviroments. Finally, we address 
the morphological dependence of the signal, which helps to 
explain some of the trends in mass-luminosity dependence. 

4.1 Luminosity dependence 

In this subsection we model the luminosity dependence of 
the signal using the g-g lensing as a function of galaxy lu- 
minosity in all 5 SDSS colors. The lenses have been divided 
into 4 luminosity bins, with the mean luminosity varying be- 
tween 0.6L* to 7 (sec figure 1). In all cases the quantity 
we analyze is AE(_R), averaged in equally spaced radial bins. 
We assume that the errors between bins are independent. 
This is a good approximation at small separations, where a 
given background galaxy contributes typically only to one 
lens galaxy (average separation between the lens galaxies is 
around 300ft _1 kpc). At larger separations this assumption 
breaks down, but since we will be mostly focusing on the 
inner regions this will not be a significant factor for the re- 
sults of the fits (although it may affect the goodness of fit, 
see below). 

The data divided in luminosity bins do not have suf- 
ficient statistical signal to allow accurate determination of 
group/cluster fraction in each bin separately. Since this frac- 
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tion is expected to decrease with luminosity (figure 2) and 
it is already low in the lowest luminosity bin we make the 
fits using two different assumptions. In the first we assume 
a — for all but the lowest luminosity bin. In the second 
we assume a from the lowest luminosity bin also applies 
to the higher luminosity bins and we use luminosity depen- 
dent Mcutoff roughly 3 times above the mass corresponding 
to that luminosity. Since a is expected to decrease with lu- 
minosity the answers we get from the two models should 
bracket the real values. We do not fit for the mean mass in 
each bin separately, but use the power law relaton between 
luminosity and mass to fit for M* and (3. 

The results from the fits for both cases are given in ta- 
ble 1 and shown in figures 11 and 12. The contour plots 
for and (3 are shown in figures 9, 10 for all bands. 
As one can see from table 1 except for the bluest band 
v! the relation between luminosity and mass is well estab- 
lished, as already shown by the SDSS team (McKay et al. 
2001). The origin for lack of correlation in v! is the mor- 
phology dependence of the signal, as discussed in more de- 
tail below. In red bands the best fitted value for (3 is sig- 
nificantly above unity, indicating that the virial mass to 
light ratio is increasing with the halo mass (or luminos- 
ity). If light is a good tracer of stellar mass (as expected 
in the reddest bands) this indicates that the star forma- 
tion efficiency is decreasing with luminosity. For example, 
in i! the fitted values are M* = (6 ± 2) x W 11 h' 1 M & , 
a = 0.2 ± 0.05 and (3 = 1.51 ± 0.15, so that the value (3 = 1 
is more than 3-<r away from the best fitted value. This gives 
M200/L ~ 30/iMq/Lq at L*, which is the lowest lumi- 
nosity that still has a detectable signal. This increases to 
A/200 / L ~ 80/iMq / Lq at 7 L+ corresponding to the bright- 
est luminosity bin. If instead M100 is used as the virial mass 
then these values should be increased by an additional 15%. 
These results however still hide significant variations be- 
tween morphological types, as discussed below. Note that 
there is a strong correlation between M* and (3. This is be- 
cause even though L* is close to the mean luminosity of 
the sample, the luminosity distribution is broad and most 
of the signal comes from the brighter end of the luminosity 
distribution. An increase in (3 gives more mass to brighter 
galaxies, which allows for a lower M* to fit the data. 

In figure 11 the best fitted model is compared to the 
data in i' band for the case where a — except in the 
lowest luminosity bin. The model provides a reasonable fit 
to the data. The reduced x 2 ls around 2.4, indicating that 
the fit is not perfect, although there is no obvious system- 
atic deviation that would indicate a clear shortfall of the 
model, suggesting instead that the observational error bars 
may be somewhat underestimated. One explanation for this 
is neglection of correlations between the bins, although this 
should only be important at large radii. Figure 13 shows 
the best fitted model assuming (3=1. For this case we find 
Mi, = 1.5 x 10 12 /i _1 Mq in i' , but the fit is less good now. 

For the case with the same value of a for all 4 bins 
we typically find M* is about 15% lower as shown in figure 
10 and in table 1, because part of the signal in the high 
luminosity bins is now accounted for by the group/cluster 
contribution. The slope (3 is reduced as well. The reduced 
X 2 reduces to 2, so this procedure does somewhat better in 
describing the data, as shown in figure 12. The value for (3 
does not change significantly (figure 10). 
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Figure 9. 68%, 95% and 99% contour plots for (3 and as 
a function of SDSS color for the case a = in all but the 
lowest luminosity bin. The results are consistent with M* ~ 
6 - 8 X 10 11 h~ 1 MQ and M/L oc L 1 / 2 except in v! where no 
clear correlation between light and mass is detected. 
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Figure 10. Same as 9 but with equal a in all luminosity bins. 



4.2 Density selection and the fraction of 
group/cluster galaxies 

In previous subsection we analyzed the luminosity depen- 
dent data to obtain the relation between the mass and lu- 
minosity around L* galaxies. In this section we combine the 
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Figure 11. Dependence of the g-g signal on luminosity using the 
4 subsamplcs in i' (figure 1) together with the best fitted model 
(solid). Dashed shows just the group and cluster contribution. We 
assume a = except for the lowest luminosity bin. 
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Figure 12. Same as figure 11 but with same a in all luminosity 
bins. 



luminosity data into a complete sample of galaxies . This 



Figure 13. Same as figure 12 but with f3 = 1. 

allows us to determine more accurately the fraction of galax- 
ies in groups and clusters, which could not be modeled for 
each luminosity bin separately. Using (5 from the previous 
subsection and the actual luminosity distribution (figure 1) 
we fit for M* and a. In i' we find M* = (6±1) x 10 11 /i" 1 M Q 
and a = 0.19 ± 0.03. These values are consistent with the 
values obtained in previous subsection for the case with a 
independent of luminosity. The resulting fits are shown in 
figure 14, where we see that the model provides a good de- 
scription of the data (reduced \ 2 ~ !)■ It is clear that both 
the central galaxy and the group/cluster components are 
needed to explain the data. 

To further test the obtained values we use the data split 
into high and low density enviroments (McKay et al. 2001). 
The separation is done using Voronoi tesselation, which as- 
signs a weight to each galaxy according to the local den- 
sity of galaxies. The galaxies are then rank ordered and di- 
vided into two equal samples. Both samples have roughly 
equal average luminosity, so under the assumption that lu- 
minosity scales with mass independent of enviroment the 
estimated M* should be approximately equal. The result- 
ing fits are shown in figure 15. For the low density sample 
the best fit gives slightly negative a = —0.07, showing that 
the low density sample is consistent with the galaxies be- 
ing only in the field. The best fitted mass in this case is 
M* = (6-7) x lO n ft _1 M0 depending on the passband. If we 
assume a = then we find M* = (4.5±1. 5) x 10 11 h' 1 M Q for 
the low density sample in i' . For the high density sample the 
values are M* = (5 - 8) x lO n ft _1 M and a = 0.45 ± 0.04. 
Assuming a = for low density sample one would expect 



* Note that averaging the luminosity bins does not give the same 
result as the average given in McKay ct al. (2001). The reason 
for the difference is the clustering correction applied to the data. 
Since some fraction of background galaxies is correlated with the 
lens they do not contribute to the lensing, so one must correct 
upwards the measured signal. This correction increases with the 



luminosity of the lens galaxy and the proper procedure is to ap- 
ply the correction to the luminosity dependent sample and then 
average the luminosity bins. Future data with better photometric 
redshift information for the background galaxies may allow one 
to eliminate this procedure by selecting only background galaxies 
that arc far from the lens. 
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Figure 14. Top: best fitted model (solid) to the complete sample 
with parameter values given in the text. Also shown are central 
galaxy contribution (long dashed) and group/cluster contribution 
(short dashed). We do not show the clustering (halo- halo) con- 
tribution, which is negligible. Bottom: contour plot of a and M*. 
We used luminosity distribution in i' for this plot. 
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Figure 16. Contour plots in c — M t plane for i' and z' . While 
the two parameters cannot be determined separately there is a 
weak dependence of on c. 



ure 16). Given that we cannot determine c independently we 
must asses its influence on M*. As seen in figure 16 there is 
some weak anti-correlation of M 4 with c. This arises because 
higher c increases the signal in the first bin, which has a sig- 
nificant contribution to the overall fit. This correlation can 
be parametrized as M* oc c~ 0,15 and since c is expected to 
be in the range between 8 and 15 this uncertainty affects M* 
by 10% at most and so does not dominate the error budget. 



a ~ 2 x 0.2 = 0.4 for the high density sample using a ~ 0.2 
obtained from the complete sample. The actual value is con- 
sistent with this, giving another confirmation that our model 
contains all the main ingredients needed to explain g-g lens- 
ing. All the reduced x 2 values are around 1, so the model 
fully explains the data. 

We can test the assumption that for a given luminosity 
the halo mass and its profile does not depend on the envi- 
roment by comparing the M* between high and low density 
samples. If the galaxies inside larger halos do not keep any 
dark matter one should see a reduction in the high density 
sample at small radii, since for 1 — a = 0.6 of all the galaxies 
there would be no individual halo attached to them. Such a 
model would give fitted M* in the high density sample about 
60% of that in the low density sample. In fact we find that 
the best fitted value for M* are comparable, which argues 
against this assumption. While the errors are still large and 
the current data does now allow us to test this hypothesis in 
more detail we conclude that there is no evidence for a signif- 
icant halo stripping or different star formation efficiency for 
galaxies within larger halos. This is consistent with the ex- 
pectations from the noninteracting cold dark matter models 
(Ghigna et al. 2000), but may not be with self-interacting 
dark matter models which predict significant mass loss in 
subhalos (Spergel & Steinhardt 2000). 

Since we have a signal as a function of radius we should 
be able in principle to place limits on the concentration pa- 
rameter c independent of the halo mass. In practice broad 
luminosity distribution and significant clustering correction 
prevents us from doing so. We have tried to do a formal fit 
to M* and c simultaneously and the resulting contour plots 
show a strong degeneracy between the two parameters (fig- 



4.3 Morphology dependence 

So far we have assumed that morphology of the lens galaxy 
does not affect the relation between luminosity and mass, 
so that we do not have to worry about their differences in 
a sample which mixes early and late types. This is a rea- 
sonable assumption in the red/infrared bands where light is 
supposed to be a reliable measure of the total stellar mass, 
assuming the morphology of a galaxy (which is believed to 
be related to the merging history) is not correlated with the 
amount of gas in a halo that cools to form stars. Observa- 
tionally this was suggested by the SDSS team (McKay et al. 
2001), showing that mass to light ratios at a fixed radius 
in the reddest bands were very similar. This analysis did 
not account for the contribution from groups and clusters 
and for observed mass- luminosity scaling. Here we assume 
mass luminosity relation based on /3 values in table 1. We 
use the sample split into late and early types that has about 
equal number of galaxies for the two types. The morpholog- 
ical classification is based on the color information and light 
concentration index and the early type sample is dominated 
by ellipticals, but also includes SO and Sa galaxies, while the 
late type sample is dominated by Sb/Sc's (see Strateva et al. 
2001 for a thorough discussion of the relation to the Hubble 
sequence). Figure 17 shows the fraction of the two types as 
a function of luminosity for the 4 bins used here. It is clear 
that at the bright end in r' , i! and z' , the sample is domi- 
nated by the early types, while in u' it is just the opposite. 
In g' the fraction is roughly independent of luminosity. Note 
that the bright end is where most of the g-g lensing signal 
is. 

We fit for M* and a in the two subsamples. The re- 
sults are given in table 2 for all 5 colors. We do not apply 
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Figure 15. Top: best fitted model to the high density sample together with the contour plot in M* — /3 plane for luminosity distribution 
in i' . Bottom: same for the low density sample. 



any internal extinction correction which should be signifi- 
cant for the late type sample in the blue bands, so for those 
the intrinsic luminosities could be even brighter. In addi- 
tion, while Petrosian magnitudes used by the SDSS should 
account for most of the light from the late types with expo- 
nential profiles they may underestimate the light from the 
early types with de Vaucouleurs profile by up to 20%. Figure 
18 shows the best fitted model together with the data. We 
find a ~ 0.07 ± 0.04 for the late type sample, showing that 
most of the late type galaxies are in the field. For the early 
type galaxies we find a ~ 0.28±0.06, so a significantly larger 
fraction of these is in groups/clusters. This is another evi- 
dence of morphology-density relation, although in this case 
we are not probing the relation as a function galaxy density 
(Dressier 1980), but rather the fraction of galaxies around L* 
that are in groups of a broad range in mass around 10 13 Mq. 
Groups of such low mass are very difficult to study observa- 
tionally and our approach is one of a few ways to study the 
galaxy membership in these halos. The errors quoted above 
are only statistical and the parameter a depends also on the 
assumed values for e, M cu toff and y g as discussed in previous 
section. We estimate the systematic error on a should be of 
order 30%. 

The values for M* differ significantly between early and 
late types in all bands. The difference is almost a factor of 
10 in u , 6 in g' , dropping to a factor of 2.5 in r' and 2 
in i' and z' . This differs somewhat from the conclusions in 
McKay et al. (2001), where it was found that M/L for the 
two types becomes equal in z' . The main reason for the dif- 
ference is that we assume M oc L fi with /3 values from table 



1 and that we apply luminosity dependent clustering correc- 
tion to the data, which boosts the brighter early type sample 
more than the fainter late type sample. However, the error 
on Mi, for the late types is large, since the signal is rather 
weak, and in i' and z' the data are still consistent with the 
assumption of equal M/L ratios for the two types. In addi- 
tion, Petrosian magnitudes may miss up to 20% of light for 
early type profiles and much less than that for late types, so 
this may bring the mass to light ratios closer by this amount. 
Theoretically one would expect M/L to be quite similar if 
i' and z' luminosity is a reliable tracer of stellar mass and if 
star formation efficiency is independent of morphology. The 
color comparisons indicate that the difference in luminosity 
between early and late types is only 20-30% from i' to K 
band (which should be closest to tracing stellar mass). We 
discuss the validity of these assumptions further below. 

Given that the M/L differs significantly between early 
and late types and that their relative contribution changes 
with luminosity we must revisit the mass-luminosity scaling 
of previous subsection. One possibility is to use g' ', where the 
two populations are in equal ratios across all luminosity bins. 
Since at a given luminosity late types are much less massive 
the dominant contribution comes from the early types, so 
the deduced /3 = 1.3 — 1.7 is the correct scaling for early 
types regardless of (3 for late types. In i' and z' the early 
types dominate at the brightest end and contribute only 
40% to the faintest bin, so here a morphological dependence 
of M/L could induce a change in /3. On the other hand, in 
these bands the late types are only a factor of 2 less massive 
for a given luminosity. This gives f3 for early types about 
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Figure 18. Top: best fitted model to the early type sample together with the contour plot in M* — a plane for luminosity distribution 
in i' . Bottom: same for the late type sample. 



0.15 lower than the value found for the overall population. 
Similar values are found also in r' and z' . We conclude that 
(5 for early types is in the range 1.4 ± 0.2 in the red bands. 

Variation of morphology with luminosity has a larger 
effect on the correlation between luminosity and mass in v! . 
Here late types dominate at the bright end and early types at 
the faint end, so the massive galaxies are actually on average 
fainter than the less massive ones. If we use f3 — 1.4 for 
early types together with the scaling of early type fraction 
with luminosity (figure 17) we obtain (3 ~ 0.5, in agreement 
with the fitted value in this band. Lack of correlation in 
this band is therefore not necessarily a consequence of v! 
light being uncorrelated with mass, but rather of a changing 
morphological fraction with luminosity coupled to a large 
difference in M/L at a given L in that band. 

The current data sample does not allow for a determi- 
nation of mass luminosity relation for late types. We have 
assumed the same value of /3 for both types here, but it 
could well be that (3 is lower for the less luminous and less 
massive late type sample, as suggested by the semi-analytic 
models of galaxy formation (Kauffmann et al. 1999, Ben- 
son et al. 2000) . If we use f3 = 1 for the late type the best 
fitted does not change much, so the difference in M/L 
between early and late type sample cannot be explained by 
this effect. 

The errors in tables 1, 2 are assumed to be gaussian 
in M^, a and (3. For M* we show the comparison in figure 
19, where the log-likelihood is plotted against the gaussian 
approximation for early and late type galaxies. We show 
both the case assuming a constant a (where the rms error 



is given by the inverse of corresponding diagonal term in 
the curvature matrix) and marginalized over a (where the 
rms error is given by the corresponding diagonal term in 
the inverse of curvature matrix). In all cases the gaussian 
approximation gives a good description of the true likelihood 
at least out to 2-sigma level. This is a consequence of the fact 
that the signal scales almost linearly with the virial mass M* 
(figure 3). 

Since for early types the group and cluster contribution 
is significant the values of M* depend somewhat on the as- 
sumed radial distribution of galaxies and dark matter, as 
well as on the halo occupation statistics. We find that if 
the group and cluster profile is shallower then fitted values 
do not change much, while if they are more concentrated 
then part of the signal at smaller radii can be attributed 
to group/cluster contribution and this lowers M*. For ex- 
ample, assuming c = 10 both for group/cluster dark matter 
and for radial galaxy distribution results in a 25% reduction 
of Mi, for early types, while the corresponding value for late 
types changes significantly less. This is probably an upper 
limit, since c = 10 is a more appropriate value for galaxy 
sized halos rather than groups and clusters, but one should 
nevertheless keep in mind this as an additional source of 
systematic error on M+. 



5 DISCUSSION 

This paper establishes a quantitative framework on how to 
analyze g-g lensing data and applies it to the SDSS sam- 
ple of 35,000 galaxies with known redshifts and luminosi- 
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Figure 17. Fraction of early types (solid) and late types (dotted) 
in each of 4 luminosity bins for all 5 colors. Also shown is the 
magnitude of L* galaxy (dashed). 
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Figure 19. Comparison between exact log-likelihood (solid) and 
its gaussian approximation assuming rms values from the non- 
linear fitting routine (dashed). Left is for late type and right for 
early type galaxies. Narrower curves are for a fixed value of a at 
the minimum, while broader curves are marginalized over a. In 
both cases, gaussian distribution gives a good description of the 
error probability distribution. 



ties. The main qualitative new feature of the model is that 
both the dark matter around the individual galaxies and 
dark matter in groups and clusters are included. Both com- 
ponents are needed to explain the observations. The two 
contributions have a different radial dependence and can be 
determined separately. This provides important constraints 
on the galaxy formation models, which must satisfy the rel- 
ative contribution from central and group/cluster compo- 
nents. We also argue that correlations between the galaxy 
and another halo that is not the one that galaxy belongs 
to can be neglected, at least on scales below l/i _1 Mpc of 
interest here. 

The main result of this paper is determination of galaxy 
virial masses and the fraction of galaxies in groups and clus- 
ters as a function of luminosity and morphology. These pro- 
vide important constraints on the galaxy formation models. 
The average virial mass M200 of an galaxy is around 
(5 — 10) x 10 11 /i _1 Mq, depending on the passband. This 
mass varies significantly with morphology and the varia- 
tion is largest in u' , where the difference between early 
and late types can be up to a factor of 10, decreasing 
to a factor of 2-3 in r' , i' and z . For example, in 1 we 
find M200 = (3.4 ± 2.1) x 10 11 h' 1 M & for late types and 
M200 = (9.3 ±2.2) x lO 11 r 1 M for early types. While the 
signal for late types is rather weak and consequently the er- 
rors are large, they are gaussian distributed and so we can 
exclude the possibility that in red bands the M* for early 
and late types are equal. 

What does this imply for the star formation efficiency 
as a function of galaxy morphology? To address this we first 
transform from i' to K band, where the luminosity is as- 
sumed to be a reliable tracer of stellar mass for a given age 
of the population (and IMF). This transformation does not 
change the results significantly, since K — i! differs by 0.2- 
0.3 magnitudes at most between the two types (Ivezic et al. 
2001). This difference is further reduced by up to 20% be- 
cause of the missed light for the de-Vaucouleurs profile. As 
a result, the difference in M/L between the two types is re- 
duced to a factor of 2. In K band the difference in M ste n ar /L 
between early and late type population can be up to a factor 
of 2 depending on the exact ages and IMF assumed, so our 
results are consistent with the assumption that star forma- 
tion efficiency for early and late type galaxies is the same. 
The errors however are still large both on the virial masses 
and stellar mass to light ratios and there could still be up to 
a factor of 2 difference. Since the current results are based 
on only 5% of final SDSS sample one should be able to place 
much better constraints on this issue in the future. 

If we adopt the late(early) type stellar mass to light 
ratio Mstoiiar/i ~ 1.5(3)h in i' we find that at L* about 
10-15% of virial mass is converted to stars. The fraction of 
baryons inside a halo should equal Clb/Q, m , which can vary 
between 10-20% for Q, m ~ 0.3±0.1 and Q b ~ 0.04±0.01. We 
see that if our halo masses are correct a significant fraction 
of baryons, up to 100%, is converted to stars in such halos. 
For early types this fraction decreases with luminosity and 
is a factor of 2 lower at 7L*. These stellar fractions can 
be accomodated in the standard models, but require a low 
matter density and/or a high baryon density, so that the 
overall baryon to dark matter fraction is sufficiently high. 

We find that M/L oc L°- 4±0 - 2 in red bands (i' and z'). 
An increase in M/L with luminosity above L* is expected 
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theoretically from semi-analytic models of galaxy formation 
(Kauffmann et al. 1999, Benson et al. 2000). Comparison 
of these results with mass tracers such as Tully-Fisher or 
fundamental plane is complicated, because these probe the 
halo at smaller radii, where the rotation velocity of dark 
matter could be different and where baryons could play a 
significant role. A detailed analysis will be presented else- 
where. The results from such analysis show that the masses 
obtained here are consistent with the CDM picture in which 
the rotation velocity drops by roughly 1.7-1.8 from the op- 
tical to the virial radius at L*. Such a drop is expected in 
CDM models, where the rotation curve peaks at a fraction of 
the virial radius and drops beyond that. However, to explain 
such a large drop one also needs a significant additional con- 
tribution to the rotation velocity from the stellar component 
and/or baryon compressed dark matter. For early types the 
scaling M oc L 15 indicates L oc w| 0i while for the same 
galaxies the Faber- Jackson relation gives L oc a 4 (Bernardi 
ct al. 2001). In this case the ratio of optical to virial circular 
velocity depends on luminosity and drops to 1.4 at 7L*. 

Comparison with other g-g lensing results is also com- 
plicated, since none of these use luminosity information, arc 
typically at a higher redshift with poorly determined dis- 
tances and may have a different morphological composition. 
The most direct comparison can be made to the study of 
early type galaxies by Wilson et al. (2001), where color in- 
formation was used to determine approximate redshifts to 
these galaxies. A direct comparison shows that their virial 
mass for early type galaxy is around 2 x 10 12 h^ 1 Mq. 
Their sample is at a higher redshift and may not be di- 
rectly comparable to ours, but the obtained value is quite 
close to our best fitted value of M* = 2.1 x 10 12 /i _1 M Q 
in g' , which is the closest to B band used there (note that 
L* = 1.1 x 10 10 /i~ 2 Lq in both samples). In their analysis 
they assume f3 = 0.5, which is inconsistent with the SDSS lu- 
minosity dependent data, so their actual value of M* could 
even be somewhat lower because j3 and M* are anticorre- 
lated. Earlier analysis of the dynamics of satellites around 
spiral galaxies gave significantly higher masses (up to a fac- 
tor of 5, Zaritsky et al. 1997) and also did not show corre- 
lation between light and mass, but a more recent analysis 
of the SDSS data shows a better agreement with g-g lensing 
results (T. McKay, private communication). At the upper 
end of the luminosity range our results can be compared 
to group and cluster velocity dispersion analysis of Girardi 
et al. (2001). Our results agree well both on the mean lumi- 
nosity of a 10 13 /i _1 Mq halo and on the scaling of mass with 
luminosity. For late type galaxies at L* our results agree well 
with the virial mass derived from theoretical models by van 
den Bosch (2001). 

There are many aspects of the analysis that could be im- 
proved upon with better statistics which will be available in 
the future. We mention some here. We have assumed a power 
law relation between mass and luminosity with a constant 
slope f3, while theoretical models suggest /3 changes with 
mass. With the current sample the signal for galaxies less 
luminous than L* is too small to be detectable and the rela- 
tion between mass and luminosity cannot be established in 
that range. As larger samples become available this regime 
will be probed as well. Fraction of galaxies in groups and 
clusters should be determined as a function of luminosity. 
Similarly, galaxy morphology dependent analysis should be 



improved by analyzing luminosity dependence of the signal 
for each morphology type and by dividing the lens galax- 
ies into several morphological classes. With better statistics 
the fraction of galaxies in groups and clusters could be deter- 
mined as a function of luminosity and not just for the overall 
sample as done here. The radial distributions of dark matter 
and galaxies in groups and clusters as well as galaxy occupa- 
tion statistics as a function of halo mass should be quantified 
better. These effects can change the values for M* and f3 by 
20-30% and remain the dominant source of systematic error. 
Better k-correction and evolution effects should be applied 
to the data to match the bright end (which receives contri- 
butions from 2 ~ 0.2) and faint end (low z) samples. This 
effect should not exceed 0.1-0.2 magnitudes in red bands and 
so should not change our conclusion on the mass luminosity 
relation significantly. Differences between halos of galaxies in 
the field and in groups and clusters could be studied in more 
detail, testing for example the collisionless CDM paradigm. 
In this paper we assumed that the profiles of galaxies within 
groups and halos are unchanged out to the truncation radius. 
In the extreme opposite case where galaxies inside groups 
and clusters do not retain any mass at all the virial masses 
of galaxies in the field will be underestimated by 1 — a. This 
is a small correction for late type galaxies, but potentially 
more important for early type galaxies. Finally, the cluster- 
ing correction, which corrects for the fraction of galaxies that 
are correlated with the lens and are not in the background, is 
luminosity dependent and should be done more accurately. 
While the procedure adopted here is reasonable, it should 
be studied in more detail. This effect may be removed for 
example if one can select background galaxies using their 
photometric redshift (photo z) information, although it is 
not clear if the faint galaxies that dominate the background 
population have sufficient signal to noise for this purpose. 
Current experience with photo z's indicates that in SDSS 
these work only down to r' < 20.5, while most of the back- 
ground galaxies used for shape information are fainter than 
that. 



As we obtain more data better statistical analysis will 
also be needed. Rather than divide the data in arbitrary 
luminosity bins, which may still be too broad for a quanti- 
tative analysis, one should parametrize the model with a few 
parameters and use maximum likelihood type of analysis to 
determine these. Photometric redshifts would also help to 
improve the statistics since one could weight the signal by 
giving optimal weight for each lens-background galaxy pair 
(this would then automatically downweight all the galaxies 
that are close to the lens). These improvements will allow 
one to make very robust statements on the relation between 
the virial mass and luminosity and on the membership in 
groups and clusters over a broad range of luminosities and 
morphological types. 
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band 


M*[lO 11 fe- 1 M ] 


13 


a 


Rm„/3 


L*[10 10 ?i- 2 LQ] 


v! 


8.05±1.23 


0.23±0.22 


0.23±0.03 


-0.04 


0.78±0.06 


g' 


10.24±1.62 


1.67±0.22 


0.20±0.03 


-0.71 


1.11±0.04 


r' 


8.96±1.59 


1.51±0.16 


0.17±0.03 


-0.83 


1.51±0.04 


i' 


8.54±1.53 


1.51±0.14 


0.17±0.03 


-0.85 


2.05±0.08 


z' 


10.01±1.66 


1.44±0.13 


0.16±0.03 


-0.83 


2.58±0.10 


u' 


6.16±1.70 


-0.52±0.29 


0.15±0.02 


0.83 


0.78±0.06 


<?' 


9.33±1.58 


1.16±0.23 


0.17±0.02 


-0.35 


1.11±0.04 


r' 


7.56±1.55 


1.34±0.17 


0.18±0.03 


-0.72 


1.51±0.04 


i' 


7.16±1.53 


1.40±0.15 


0.17±0.03 


-0.79 


2.05±0.08 


z< 


8.08±1.62 


1.35±0.14 


0.17±0.03 


-0.75 


2.58±0.10 



Table 1. Best fitted parameters to the SDSS data in five colors assuming a = for all but the faintest bin. The corresponding values for 
Mi, using equal a for all bins are 10-20% lower, while the other parameters do not change much. Also given is the correlation coefficient 
between M, and and the value of L*. 



band 


type 


M t [W 11 h- 1 M Q ] 


a 


Li,[W 10 h- 2 L Q ] 


u' 


E 


25.55±5.78 


0.31±0.07 


0.78±0.06 


u' 


S 


2.52±1.58 


0.09±0.04 


0.78±0.06 


g' 


E 


21.94±5.25 


0.29±0.07 


1.11±0.04 


g' 


S 


3.18±2.02 


0.08±0.05 


1.11±0.04 


r' 


E 


10.73±2.53 


0.28±0.08 


1.51±0.04 


r< 


S 


3.26±2.08 


0.08±0.05 


1.51±0.04 


i' 


E 


9.32±2.26 


0.28±0.08 


2.05±0.08 


i' 


S 


3.44±2.14 


0.08±0.05 


2.05±0.08 


z' 


E 


10.06±2.38 


0.28±0.08 


2.58±0.10 


i 


S 


4.08±2.63 


0.08±0.05 


2.58±0.10 



Table 2. Best fitted parameters for early type (E) and late type (S) galaxies as a function of color. Early types have higher M* and 
group/cluster fraction than late types. 
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